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An  unexpected  mechanism  for  fast  oxidation  of  A1  nanoparticles  covered  by  a  thin  oxide  shell  (OS) 
is  proposed.  The  volume  change  due  to  melting  of  A1  induces  pressures  of  0.1-4  GPa  and  causes 
spallation  of  the  OS.  A  subsequent  unloading  wave  creates  high  tensile  pressures  resulting  in 
dispersion  of  liquid  A1  clusters,  oxidation  of  which  is  not  limited  by  diffusion  (in  contrast  to 
traditional  mechanisms).  Physical  parameters  controlling  this  process  are  determined.  Methods  to 
promote  this  melt  dispersion  mechanism,  and  consequently,  improve  efficiency  of  energetic 
nanothermites  are  discussed.  ©  2006  American  Institute  of  Physics.  [DOI:  10.1063/1.2335362] 


Thermites  are  mixtures  of  metal  fuel  (e.g.,  Al,  Mg,  and 
B)  and  oxidizers  (e.g.,  Fe203,  Mo03,  and  CuO)  which  react 
to  produce  high  temperatures  with  slow  oxidation  reactions. 
Since  fuel  particles  are  covered  by  a  thin  OS,  oxidation  is 
controlled  by  diffusion  of  molecules  through  a  growing  OS. 
However,  when  the  initial  external  radius  R  of  the  fuel  par¬ 
ticle  is  reduced  to  10-60  nm  (in  contrast  to  traditional 
1  — 100  /xm  size),  reaction  rates  drastically  increase  and 
flame  propagation  rates  C  reach  0.9-1  km/s,  1-3  in  contrast 
to  cm/s  for  traditional  thermites.  Rise  time  for  temperature 
and  pressure  at  the  reaction  front  is  tf—  10  /zs  which  charac¬ 
terizes  the  reaction  time.  For  self-diffusion  coefficients  of 
oxygen  and  Al  in  a-alumina  at  800-950  °C,  D- 10~19  and 
10-18  cm2/s,4  the  diffusion  length  ld=2^Dtf=(2-6) 
X  10-5  nm,  which  is  105  times  smaller  than  the  OS  thickness. 
Also,  C  and  ignition  time  delay  tig  are  independent  of  par¬ 
ticle  size  for  R<  40-50  nm,  while  for  diffusion  controlled 
oxidation  C—l/R  (Ref.  2)  and  tig~Rn.5  Thus,  the  only 
known  oxidation  mechanism  via  diffusion  does  not  apply  for 
nanoparticle  combustion;  finding  such  a  mechanism  for 
nanosize  thermites  is  one  of  the  most  important  problems  in 
combustion  physics.  In  this  letter,  we  suggest  and  justify 
theoretically  and  experimentally  an  unexpected  oxidation 
mechanism  during  fast  heating  which  explains  numerous 
puzzles  observed  experimentally  in  flame  propagation  of 
A1+Mo03  (Refs.  1,  2,  and  6)  and  Al+Fe203  (Ref.  3)  nano¬ 
thermites.  Fast  heating  creates  huge  thermal  stresses  due  to  a 
difference  in  thermal  expansion  coefficients  a  of  Al  and 
A1203  and  volumetric  strain  3sm  during  the  Al  melting.  Due 
to  the  small  size  of  the  amorphous  OS,  £=1-4  nm,  it  is 
defect-free  and  its  strength  approaches  the  theoretical 
strength.  For  small  particles  (small  M=  R/ S)  the  OS  frac¬ 
tures  after  melting  of  an  Al  particle  (in  contrast  to  slow  heat¬ 
ing  or  micrometer  size  particles).  Pressure  in  the  melt  is  very 
high  (0.1-4  GPa).  Dynamic  spallation  of  OS  results  in  com¬ 
plete  exposure  of  the  liquid  Al  droplet  which  leads  to  an 
unloading  wave  within  the  liquid  Al  droplet  and  creates  high 
tensile  pressure.  The  tensile  pressure  disperses  the  Al  droplet 
into  small  clusters  that  are  ejected  at  a  high  velocity.  The 
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oxidation  of  these  clusters  is  not  limited  by  diffusion.  The 
suggested  melt  dispersion  mechanism  (MDM)  is  much  faster 
than  diffusion  and  is  independent  of  the  size  of  Al  particle 
(below  some  critical  size);  fast  reaction  provides  fast  flame 
propagation  and  fast  heating  rates  for  the  next  particles.  The 
necessary  conditions  for  operation  of  this  mechanism  and 
controlling  physical  parameters  are  found.  By  altering  con¬ 
trolling  physical  parameters  we  suggest  ways  to  promote  the 
MDM,  and  thereby  improve  the  efficiency  of  energetic  nano¬ 
thermites.  Increasing  the  temperature  T0  at  which  initial  OS 
is  formed,  creating  initial  porosity  in  the  Al,  and  alloying  Al 
to  decrease  the  cavitation  pressure  are  ways  to  induce  the 
MDM.  Theoretical  results  are  consistent  with  observations. 
A  similar  mechanism  is  expected  for  nitridation  and  fluoriza- 
tion  of  Al.  The  MDM  may  also  be  tailored  for  Ti  and  Mg 
fuels. 

Let  us  consider  internal  stresses  during  the  heating  of 
spherical  Al  particle  with  radius  R  covered  by  a  thin  alumina 
OS  with  external  radius  R.  The  surface  temperature  of  the 
particle  varies  from  Tr=  300  K  to  the  final  Tf—  1000  K  (just 
above  melting  temperature  of  Al,  Tm= 933.67  K)  during  a 
time  tf.  It  is  easy  to  show  that  for  R  <  9  jam  the  temperature 
in  the  Al  sphere  is  homogeneous.  Stresses  for  a  large  Al 
particle  covered  by  OS  were  found  in  Ref.  7;  we  generalize 
these  results  for  small  particles  and  apply  the  proper  strength 
criterion.  Internal  stresses  arise  during  heating  due  to  a  dif¬ 
ference  in  inelastic  strains  sl2-s\  (subscript  1  is  for  Al  and  2 
is  for  A1203).  Surface  tensions  at  the  A1-A1203  interface,  T  b 
and  between  alumina  and  gas,  T 2,  cause  jumps  in  pressures 
of  2 Tl/R  (at  r=R )  and  2 T2/R  (at  r=R ).  Pressure  in  the  Al 
sphere,  p ,  and  maximum  hoop  stresses  in  OS  at  r=R ,  ah ,  can 
be  found  from  the  elasticity  theory:8 


12(m3  -  l)(e2  -  s[)G2KiK2  2Kl(4G2  +  3m3K2)Tl 

H  +  RH 

2m2Kl(4G2  +  3K2)T2 
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TABLE  I.  Material  parameters  at  T=Tm. 


K\ 

K " 

k2 

g2 

* 

“2 

r 

P\ 

YYl 

Pi 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

(105  K->) 

(105  K->) 

(105  K"1) 

sm 

(GPa  nm) 

(kg/m3) 

(kg/m3) 

71. la 

41.3b 

234.8C 

149.5C 

3.032a 

4.268d 

0.778c 

0.02e 

1.05f 

2530g 

2380d 

“See  Ref.  13. 
bSee  Ref.  14. 
“See  Ref.  12. 
dSee  Ref.  16. 


“See  Ref.  7. 
fSee  Ref.  17. 
gSee  Ref.  15. 
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2 m2(-  2 G2Kx  +  3(2 G2  +  Kl)K2)T2 

RH 


RH 


(2) 


e[  =  ~  (a\(Tm  -  T0)  +  (1  -f)a\{T- Tm)  Tm) 

+  /8-),  (3) 

s2  =  -  a2 (T-  To). 

Here  K  and  G  are  the  bulk  and  shear  moduli,  Kx  =fK™+{\ 
-f)K\,  m=R/R=  1  +  S/R ,  /  is  the  volume  fraction  of  the  melt 
in  Al  particle,  subscripts  s  and  m  are  for  the  solid  and  melt 
phases,  and  H=3m3KlK2+4G2(Kl  +  (m3-  l)K2).  The  crite¬ 
rion  for  brittle  fracture  of  OS  (maximum  tensile  stress 
reaches  the  ultimate  tensile  strength  au )  is  -crh=cru. 

All  material  parameters  are  given  in  Table  I  (T1  =  r2 
=  T).  For  amorphous  OS  of  £=1.5-8  nm,  defects  are  ne¬ 
glected,  thus  the  ultimate  strength  is  equal  to  theoretical 
strength,  ath.  For  crystalline  materials,  ath=E/20-E/\0  but 
we  conservatively  take  au=E2/3 0=11.33  GPa  and  consider 
much  smaller  cru  as  well.  Also,  because  of  lack  of  stress 
concentrators  and  fast  loading,  fracture  and  spallation  occur 
throughout  the  entire  OS. 

Using  Eq.  (2)  and  applying  the  condition  -crh=cru,  the 
value  of  volume  fraction  of  melt,  £,T0),  necessary 

for  the  fracture  of  the  OS  is  obtained.  The  functions  fj(M) 
for  r0=600  K  and  several  values  of  cru  and  S=  3  are  shown 
in  Fig.  1(a).  They  do  not  change  significantly  for  1.5  <S 
<  8  nm.  In  small  particles  the  fracture  criterion  is  not  ful¬ 
filled  even  at  complete  melting.  Then  OS  spallation  will  oc¬ 
cur  at  T>  Tm,  due  to  growth  of  thermal  stresses  and  reduc¬ 
tion  in  cru  with  the  temperature  increase.  In  such  a  case,  an 
unloading  wave  in  Al  particle  after  OS  spallation  produces 
tensile  pressure  and  disperses  the  liquid  Al  droplet  (see  be¬ 
low). 

Relationships  Po(M)  for  pressure  in  Al  particle  at  the 
instant  of  OS  fracture  are  shown  in  Fig.  1(b).  For  cru 
=  9.0  GPa  and  M=19.0,  pressure  at  the  complete  melting  is 
1.0  GPa;  for  smaller  cru ,  complete  melting  occurs  at  smaller 
M  and  pressure  is  even  larger.  For  comparison,  for  T0 
=  300  K  we  obtain  complete  melting  at  M=10.8  and  p 
=  1.8  GPa;  for  T0=900  we  have  M- 80  and  p  =  0.24  GPa. 
Thus,  T0  is  an  important  parameter  controlling  which  MDM 
can  be  produced  even  for  large  particles. 

Thus  at  fast  heating,  the  major  part  of  the  Al  nanopar¬ 
ticle  melts  before  fracture  of  OS,  producing  high  pressure. 
Then  unloading  wave  after  OS  spallation  disperses  the  liquid 
particle.  For  smaller  particles,  a  larger  portion  of  the  Al  melts 
and  the  pressure  in  the  Al  is  higher  (Fig.  1).  The  MDM 


should  not  occur  during  slow  heating9  and  for  large  particles. 
For  cru= 9.0  GPa  and  R/ S=  2000,  we  obtain  p0<9.0MPa 
which  may  not  produce  sufficient  tensile  stress  in  the  unload¬ 
ing  wave  to  disperse  solid  Al  particles.  For  T=Tm,  cru 
=  9.0  GPa,  and  T0  =  300  K,  we  have/f=0.21  only. 

Let  a  sphere  of  radius  R  be  initially  in  equilibrium  under 
applied  external  pressure  p0.  We  found  analytically  pressure 
and  velocity  evolutions  after  reduction  of  the  pressure  at  the 
boundary  r=R  from  p0  to  the  final  value  Pf  during  the  time 
of  OS  spallation,  ts,  after  which  it  remains  constant.  We  de¬ 
fine  an  acoustic  time  tp=R/c  during  which  boundary  regime 
propagates  to  the  center  of  the  sphere,  where  c=  p™ 
=  4166  m/s  is  the  sound  velocity  in  Al.  Thus,  for  R 
=  41.66  nm  we  have  tp=\0  ps.  We  consider  unloading  and 
the  first  reflected  wave  only.  For  fast  heating  of  an  Al  nano¬ 
particle  with  a  thin  OS  when  the  OS  fractures  at  reaching  its 
theoretical  strength,  ts<tp\  we  take  ts=0.2tp.  For  R  =  40  nm, 
we  take  p{)-  1  GPa  and  pf=2T/R  —  0.05  GPa. 

The  results  shown  in  Fig.  2  are  typical  for  the  regime 
with  ts<tp  and  Pf^Po-  For  ts<t<tp,  a  region  with  tensile 
pressure  appears  and  propagates  to  the  center  of  the  sphere; 
the  magnitude  of  tensile  pressure  grows  drastically  (Fig.  2). 
Maximum  values  of  tensile  pressure  and  velocity  in  the  cen¬ 
tral  part  of  the  sphere  reach  |/?|  =  3.75  GPa  and  u  =  240  m/s, 
respectively. 

Fracture  of  liquid  (cavitation)  criterion  is  p<pc< 0, 
where  pc  is  the  critical  tensile  pressure;  we  can  estimate  pc 


FIG.  1.  Volume  fraction  of  melt  necessary  for  fracture  of  OS  ff  (a)  and 
pressure  p0  in  Al  particle  at  the  fracture  of  OS  (b)  as  a  function  of  M 
=R/  S  for  various  values  of  ultimate  strength  cru  (shown  near  the  curves); 
<5=3  nm,  r0=600  K.  The  upper  line,  common  for  all  curves,  corresponds  to 
complete  melting.  Diamonds  represent  dimensionless  flame  velocity  C 
=  C/(  1  km/s)  vs  M  from  experiments  (Ref.  2).  Triangles  are  the  relative 


^  v  ,  ignition  time  ?ig/ ( 1 00  ms)  from  experiments  (Ref.  5). 
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FIG.  2.  Pressure  (GPa)  distribution  in  Al  sphere  after  spallation  of  OS  at 
different  instances  of  propagation  of  the  unloading  wave  (shown  near  the 
curves  as  a  fraction  of  tp).  The  upper  line  (t=0,p=l  GPa)  and  the  lowest 
curve  are  common  for  all  curves  for  t^tp,  i.e.,  at  each  time  the  pressure 
profile  consists  of  straight  horizontal  line  ( 1  GPa)  in  the  central  part,  the  line 
connecting  to  the  lowest  curve  and  part  of  the  lowest  curve. 

=  -(0.01-1)  GPa.  For  pc=- 0.1  GPa,  the  cavitation  criterion 
is  satisfied  at  radius  7y=  0.857?  (Fig.  2).  Cavitation  results  in 
separation  of  a  thin  liquid  outer  shell  from  the  main  droplet. 
The  thickness  of  this  shell  during  flight  will  get  smaller 
down  to  atomic  size  then  it  will  break  into  small  sectors. 
Since  tensile  pressure  without  fracture  exceeds  significantly 
\pc\  (Fig.  2),  the  thickness  of  each  subsequent  separated  liq¬ 
uid  shell  is  of  order  of  several  atomic  sizes.  Al  clusters  may 
react  with  oxygen  and  nitrogen  in  the  air,  gaseous  Mo03 
which  starts  to  sublimate  at  973  K,  with  Teflon  which  de¬ 
composes  at  825  K,  or  solid  oxidizers. 

For  the  same  conditions  but  for  ts=tp,  there  is  no  tensile 
pressure  in  Al  for  t<tp.  For  t>tp  cavitation  is  in  principle 
possible  but  not  very  probable.  For  ts  =  2tp,  tensile  pressure 
does  not  appear  in  the  unloading  and  first  reflected  waves, 
i.e.,  dispersion  cannot  occur. 

If,  e.g.,  cru= 9  GPa,  then  for  19  the  entire  particle 
melts  before  OS  fracture  [Fig.  1(a)]  (pressure  p0=  1.0  GPa 
for  £=3)  and  disperses  providing  maximum  possible  mixing, 
oxidation  rate,  energy  release,  and  flame  front  velocity.  Ac¬ 
cording  to  the  MDM,  there  should  not  be  an  essential  differ¬ 
ence  in  the  mixing  of  components  and  oxidation  rate  of  par¬ 
ticles  with  R/ 19.  In  our  experiments,2  flame  velocity  is 
approximately  the  same  (950  m/s)  for  the  Al  particles  with 
an  averaged  R  =  20  and  38  nm  and  reduces  for  particles  with 
R  =  53  and  58.5nm  (S—  2  nm)  [Fig.  1(a)].  In  our  recent 
experiments,5  laser  induced  ignition  times  are  in  the  range  of 
12.4-25.5  ms  for  7?=6.7-99  nm  and  grows  to  6.039  s  for 
R  =  20  /z m  [Fig.  1(a)].  These  results  support  our  theory. 

Damage  of  the  OS  should  suppress  the  MDM  but  pro¬ 
mote  traditional  diffusion  oxidation  mechanism.  Pressing 
thermite  pellets  in  a  die  potentially  damages  the  OS.  That  is 
why  C  is  expected  to  decrease  with  sample  density  for  nano¬ 
thermites  and  grow  for  micro  thermites.  That  is  exactly  what 
was  observed  in  our  experiments.6  On  the  other  hand,  vibra¬ 
tional  compaction  without  pressure  should  not  damage  OS 
and  suppress  the  MDM.  Indeed,  vibrational  change  in  den¬ 
sity  in  the  range  of  5%-10%  of  TMD  did  not  change  C.2 

Consider  an  ellipsoidal  penny-shaped  particle  with  semi¬ 
axes  b  and  R  >  b  that  mimics  a  nanoflake.  The  hoop  stress  in 
the  circular  section  of  radius  R  is  a^ax  —  -pR/ (2  S),  while  the 
hoop  stress  in  the  elliptic  section,  o^m  —  -pirbl (AS),  is  much 
smaller.  Simultaneous  spallation  of  the  entire  OS  cannot  oc¬ 
cur.  Fracture  occurs  along  the  circular  section  of  radius  R 
when  a^ax=-crM,  i.e.,  for  p  =  2auS/R  [like  for  large  spherical 
particles].  For  nanoflakes,  Z?=10-100nm,  while  R 
=  10-100  /zm.10  Thus,  fracture  of  a  flake  occurs  at  a  similar 
pressure  as  a  micro  particle.  Thus,  according  to  the  MDM, 


nanoflakes  would  tend  to  react  like  micro  particles,  while 
traditional  arguments  (higher  reactivity  and  large  surface  to 
volume  ratio)  suggest  similar  reactivity  for  nanoflake  and 
nanospherical  particles.  In  experiments,10  nanoflakes  did  not 
produce  higher  flame  rate  than  micro  spherical  particles. 

The  main  conditions  for  the  dispersion  of  the  spherical 
Al  particle  are  that  it  completely  melts  before  OS  spallation, 
that  OS  fracture  time  is  smaller  than  wave  propagation  time 
(ts  <  tp ),  and  that  pressure  in  the  particle  before  OS  spallation 
is  approximately  equal  to  the  cavitation  pressure.  These  con¬ 
ditions  can  be  fulfilled  for  a  thin  amorphous  defect-free  or 
nanocrystalline  OS,  and  only  during  fast  heating.  We  can 
predict  several  ways  to  increase  particle  radius  for  which  the 
MDM  is  operative. 

(1)  Our  analysis  predicts  that  increasing  T0  may  extend  the 
MDM  for  large  particles.  If  we  increase  T0  up  to  900  K, 
then  for  au  =  9.0  GPa  almost  the  whole  particle  melts  for 
any  particle  radius.  If  we  assume  that  Po=20  MPa  is 
sufficient  for  fracture  and  dispersion  of  Al  particle,  then 
for  au  =  9.0  GPa  and  S=  8  nm,  we  obtain  R  =  7.2  /zin . 

(2)  To  create  the  condition  for  fast  heating  of  large  particles, 
they  have  to  be  mixed  with  nanosized  particles  (7? 
<40-100  nm)  which  have  much  higher  probability  to 
be  dispersed  in  a  liquid  state.  Indeed,  a  mixture  of  70% 
of  Al  particles  with  R  =  38  nm  (R/ S=  8.3)  and  with  30% 
particles  with  R  =  2  and  10  /zm  provided  practically  the 
same  flame  velocity  as  100%  of  Al  nanoparticles.11 

(3)  An  alternative  way  to  decrease  internal  stresses  before 
melting  and  to  ensure  complete  melting  before  OS  frac¬ 
ture  is  to  introduce  initial  porosity  in  the  Al  particles. 

(4)  Alloying  the  fuel  particle  may  reduce  its  local  surface 
energy  and  cavitation  pressure. 
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Fast  oxidation  of  A1  nanoparticles: 
diffusion  models  do  not  describe  this  process 


-  Rise  time  for  T  and  P  T~10  ps 

•  characterizes  the  reaction  time 

•  Diffusion  coefficient  of  Oxygen  &  A1  in  a- 
alumina  at  800  °C  is  D  ~10'19  cm2/s 

-  Diffusion  length  ld  =  2«jL)T  =  (2-6)1  O'5  nm 
105  times  smaller  than  oxide  shell  thickness!! 

•  Diffusion  time  would  need  to  be  increased  by 
a  factor  of  10 10  to  obtain  a  reasonable 
diffusion  length  of  2-6nm. 
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-  Burn  rate  is  independent  of 

particle  size  for  dA1  <100  nm,  while  for 
diffusion  controlled  oxidation  ~l/r 

-  Ignition  delay  time  independent  of  particle 
size  for  dA1  <100  nm,  while  for  diffusion 
controlled  oxidation  t~rn 

-  Nanoflakes  react  like  micron  size  particles 

-  Previous  diffusion  models  do  not  explain 
experimental  observations ! 


A  TEM  of  an  Al  particle 


New  Mechanism  for  Fast  Reactions 

•Due  to  small  size  of  the 
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Atomic  size  molten 
aluminum  clusters 
disperse  from  an 
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amorphous  shell  and  very  fast 
m  strain  and  stress  growth 
(~106GPa)  strength  approaches 
theoretical  strength. 

•Pressure  in  the  melt  is  very  high 
(0.1-4  GPa) 

•Dynamic  spallation  of  shell 
results  in  complete  exposure  of 
the  liquid  Al  droplet  and  creates 
a  tensile  pressure  (0.1-4  GPa) 

•The  tensile  pressure  in  an 
unloading  wave  disperses  the  Al 
droplet  into  small  clusters  which 
fly  with  high  velocity 

•Oxidation  is  not  limited  by 
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Exploit  the  melt  dispersion  mechanism 


•  Suggest  methods  that  will  improve  the  performance  of  nano 
and  micron  size  composites  by  exploiting  this  new 
mechanism  for  fast  oxidation. 

•  Reduction  of  sensitivity  to  the  ignition  in  combination  with 
high  performance  (flame  rate)  is  the  main  concern. 

Controlling  parameters 

(1)  the  temperature  at  which  the  oxide  shell  is  formed  determines 

initial  stresses  in  the  particle  and  core; 

(2)  the  porosity  of  the  A1  core; 

(3)  ratio  of  oxide  thickness  to  particle  diameter; 

(4)  defects  of  the  oxide  shell; 

(5)  particle  size  distribution,  including  combination  of  nano-  and 

micron  size  particles. 
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